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CD 

polar  coordinates  centered  at  the  center  of  the  tube 
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radius  to  the  tip  of  a  crack,  rc  =  1  +  c 
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wall  thickness  of  the  cylinder,  t  =  b  -  1 

T,T0,Tp 

temperature  at  r,  r  =  a  and  r  =  p,  respectively 
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linear  thermal  expansion  coefficient 
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percentage  of  overstrain,  e  =  (p~l)/t,  0  <  e  <  1 
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Poisson’s  ratio 
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radius  of  elastic-plastic  interface  during  pressurization 

°o 

uniaxial  yield  stress  of  the  tube  material 
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CD 

normal  stress  in  the  radial  and  tangential  direction, 
respectively 
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INTRODUCTION 


In  a  previous  paper, 1  stress  intensity  factors  were  obtained,  using 
12-node  quadrilateral,  isoparametric  elements,  for  a  uniform  array  of  equal 
depth  radial  cracks  originating  at  the  internal  boundary  of  a  pressurized 
thick-wall  cylinder.  To  increase  the  maximum  pressure  a  cylinder  can  contain, 
it  is  a  common  practice  to  produce  a  favorable  residual  stress  in  the  cylinder 
by  an  autofrettage  process.  It  is  important  to  find  the  effect  of  residual 
stresses  on  the  stress  intensity  factor  for  a  cylinder  with  multiple  cracks. 

Using  the  concept  of  thermal  simulation, 2  the  autofrettage  residual 
stresses  are  simulated  by  active  thermal  loads.  It  is  shown  that  the  stress 
intensity  factor  for  multiple  radial  cracks  in  a  tube  with  residual  stresses 
can  be  computed  by  the  same  finite  element  method. 

Slight  changes  in  geometrical  configurations  or  loading  conditions 
require  new  computations.  To  obviate  this  problem,  which  is  a  shortcoming  of 
the  finite  element  method,  load  relief^  and  weight  function  methods^  are 
examined.  For  a  small  number  of  radial  cracks,  the  method  of  load  relief 
enables  us  to  estimate  the  stress  intensity  for  N  other  than  two  fairly  accu- 


*Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  For  a  Circular  Ring 
With  Uniform  Array  Radial  Cracks  Using  Isoparametric  Singular  Elements,"  ASTM 
STP-677 ,  1979,  pp.  685-699. 

^Hussain,  M.  A.,  Pu.  S.  L.,  Vasilakis,  J.  D.,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 

Baratta,  F.  I.,  "Stress  Intensity  Factors  For  Internal  Multiple  Cracks  in 
Thick-Walled  Cylinders  Stressed  by  Internal  Pressure  Using  Load  Relief 
Factors,"  Engineering  Fracture  Mechanics,  Vol.  10,  1978,  pp.  691-697. 
^Bueckner,  H.  F.,  "A  Novel  Principle  For  the  Computation  of  Stress  Intensity 
Factors,"  Z.  Agnew.  Math.  Mech. ,  Vol.  50,  1970,  pp.  529-546. 

JRice,  J.  R. ,  "Some  Remarks  on  Elastic  Crack-Tip  Stress  Fields,"  Int.  Journal 
of  Solids  and  Structures,  Vol.  8,  1972,  pp.  751-758. 
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rately  by  making  use  of  the  finite  element  result  for  N  =  2.  When  N  or  crack 
depth  is  large  it  is  shown  in  this  report  that  the  load  relief  factor  is  not 
reliable  because  it  varies  with  the  nature  of  the  load. 

For  a  given  geometrical  configuration  it  is  useful  to  use  the  weight 
function  method.  In  the  present  approach  the  restrictive  assumption  that  the 
crack  opening  displacement  is  a  conic  section® is  circumvented.  There 
are  only  three  types  of  hoop  stress  namely:  constant,  1/r2  and  log  (r) ,  in  an 
uncracked  cylinder  subjected  to  internal  pressure,  uniform  tension  on  the 
outer  boundary,  and  the  autofrettage  residual  stress.  This  eliminates  the 
need  to  assume  the  crack  face  pressure  as  a  simple  polynomial.®  It  enables 
us  to  obtain  stress  intensity  factors  for  each  type  of  crack  face  pressure 
from  three  linear  algebraic  equations  using  three  finite  element  results  for 
the  given  geometry.  The  stress  intensity  factor  can  be  readily  calculated  for 
any  combination  of  internal  pressure  and  any  degree  of  autofrettage. 

When  the  radial  cracks  progress  beyond  the  elastic-plastic  interface 
produced  during  the  autofrettage  overstrain,  the  algebraic  equation  for  stress 
intensity  factor  breaks  down  because  the  crack  face  pressure  cannot  be 
represented  by  a  simple  expression.  Modifications  are  derived  for  such  cases 
based  on  the  crack  opening  displacement  near  a  crack  tip  being  parabolic. 9 


®Orange,  T.  W. ,  "Crack  Shapes  and  Stress  Intensity  Factors  for  Edge-Cracked 
Specimens,"  ASTM  STP-513,  1972,  pp.  71-78. 

^Grandt,  A.  F.,  "Stress  Intensity  Factors  For  Some  Through-Cracked  Fastener 
Holes,"  Int.  Journal  of  Fracture,  Vol.  11,  1975,  pp.  283-294. 

®Grandt,  A.  F.,  "Stress  Intensity  Factors  For  Cracked  Holes  and  Rings  Loaded 
With  Polynomial  Crack  Face  Pressure  Distributions,"  Int.  Journal  of  Fracture, 
Vol.  14,  1978,  pp.  R221-R229 . 

^Paris,  P.  C.  and  Sih,  G.  C.,  "Stress  Analysis  of  Cracks,"  ASTM  STP-381,  1965, 
pp.  39-81. 
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Extensive  numerical  results  are  presented  for  multiply  cracked  cylinders 
having  an  external  diameter  twice  that  of  the  internal  diameter.  Results 
agree  with  those  obtained  by  other  methods  for  the  non-autofrettaged^  and 
fully  autof rettagedll  cases.  In  addition,  results  are  given  for  the 
partially  autof rettaged  cases. 


RESIDUAL  STRESS  AND  THERMAL  SIMULATION 


For  the  case  of  plane  strain,  the  stress  distribution  of  a  partially 


autof rettaged  tube,  using  the  von  Mises'  yield  criterion  for  the  incompress¬ 


ible  material,  is  given  by^2 


(  0 o  r  p2  11 

—  {(2  log  -  -  1  +  -t)  -  PiC— 5  -  -o)}  1  <  r  <  p  (1) 

/T  p  bz  bz  rz 

or(r)  =  < 


°o  ,  11 

—  (p2-Pi) (--  -  --)  P  <  r  <  b  (2) 

/3  b2  r2 


°o  r  p2 
—  2  log  -  +  1  +  -- 
/3  P  b2 


°e(r)  =  { 


°o  11 

—  (p2-Pi)(-r  + 

/3  b2  r2 


1  1  i 

-  Pi(--5  +  -^)  } 


1  <  r  <  p 


p  <  r  <  b 


(3) 

(4) 


l^Tracy,  P.  G.,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer 
and  Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol. 
11,  1979,  pp.  291-300. 

^Parker,  A.  P.  and  Andrasic,  C.  P.,  "Stress  Intensity  Prediction  For  a 

Multiply-Cracked,  Pressurized  Gun  Tube  With  Residual  and  Thermal  Stresses," 
Presented  at  Solid  Mechanics  Symposium,  Cape  Cod,  MA,  1980. 
l^Hill,  R. ,  The  Mathematical  Theory  of  Plasticity,  Oxford  at  the  Clarendon 
Press,  1950. 
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where 


b2  p‘ 

Pi  =  Pi(p)  =  (1  -  —  +  2  log  p) 

b2-l  b2 


(5) 


If  the  same  hollow  cylinder  is  subjected  to  a  thermal  load 

(T0-Tp) 


To  - 


T(r)  = 


log  p 


log  r 


1  <  r  <  p 


p  <  r  <  b 


the  thermal  stresses  are  given  by^ 
f  Ea(TQ-Tp) 


or(r)  =  \ 


r  r  P  1  1  i 

{(2  log  -  -  1  +  --)  -  PX(--  -  —J]  1  <  r  <  p 

2(l-v)log  p  p  b2  b2  r2 


Ea(T0-Tp)  i  i 

- -  (p2_pl)(  _  } 

2(l-v)log  p  b 2  r2 


Ea(X0-Tp) 


J0(r)  =  { 


r  r  P  A  1  i 

{2  log  -  +  1  +  --  -  Pi(--  +  ~r) }  1  <  r  <  p 

2(l-v)log  p  p  b2  b2  r2 


Ea(T0  Tp) 

"7 - 7 -  (P2-Pl)(~^  +  -r) 

2 ( l-v)log  p  b2  r2 


(6) 


(7) 


P  <  r  <  b  (8) 


(9) 


P  <  r  <  b  (10) 


The  thermal  stresses  and  the  autofrettage  residual  stresses  become  equivalent 
if  the  temperature  gradient  of  the  thermal  load  and  the  uniaxial  yield  stress 
of  the  cylinder  material  have  the  following  relation 

Ea(T0-Tp)  2a0 


2(l-v)log  p  /3 

where  T0  or  Tp  may  be  assigned  arbitrarily. 


(ID 


2Hussain,  M.  A.,  Pu,  S.  L.,  Vasilakis,  J.  D.,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 
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In  many  instances  there  is  a  redistribution  of  residual  stresses  due  to 
changes  of  geometrical  configurations  such  as  the  presence  of  keyways,  holes, 
riflings,  and  cracks.  The  stress  redistribution  may  be  difficult  to  find.  In 
case  of  cracks,  it  may  be  difficult  to  obtain  stress  intensity  factors.  For 
simple  problems,  the  method  of  superposition  may  be  used  to  compute  the  stress 
redistribution  as  well  as  the  stress  intensity  factors  involving  cracks. 
Several  examples  are  given  in  reference  13.  An  alternative,  the  thermal 
simulation  method  is  proposed  in  reference  2  replacing  the  residual  stress  by 
an  active  thermal  load.  The  method  of  thermal  simulation  shows  in  reference 
13  that  the  thermal  stress  redistribution  is  equivalent  to  the  residual  stress 
redistribution.  For  more  complicated  problems,  the  thermal  simulation  method 
has  the  advantage  over  the  method  of  superposition  by  eliminating  the 
computation  of  residual  stresses  acting  on  the  surface  of  the  crack  which 
would  be  present  in  the  uncracked  body  under  external  loads.  In  this  report 
the  thermal  simulation  method  is  used  to  compute  stress  intensity  factors  for 
multiple,  radial  cracks  in  a  fully  autofrettaged  tube. 


FINITE  ELEMENT  METHOD 

The  finite  element  technique  has  become  an  important  numerical  method  for 
practical  problems  in  structural  mechanics  because  of  its  ability  to  treat 
very  general  geometrical  configurations  and  loading  conditions.  The  trend  is 


2Hussain,  M.  A.,  Pu.  S.  L.,  Vasilakis,  J.  D.,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 
lJPu,  S.  L.  and  Hussain,  M.  A.,  "Residual  Stress  Redistribution  Caused  by 
Notches  and  Cracks  in  a  Partially  Autofrettaged  Tube,  Technical  Report 
ARLCB-TR-81005 ,  Benet  Weapons  Laboratory,  LCWSL,  ARRADCOM,  US  Army,  January 
1981. 
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to  use  high  order  elements  requiring  a  small  number  of  elements  for  high 
degree  of  accuracy.  The  plane  problem  of  a  uniform  array  of  equal  depth 
radial  cracks  emanating  from  the  bore  of  a  pressurized,  non-autofrettaged  tube 
has  been  solved  using  12-node  quadrilateral,  isoparametric  elements. *  The 
collapsed  singular  elements  developed  by  Pu  and  Hussain^  are  used  around  the 
crack  tip.  The  finite  element  results  of  stress  intensity  factors  agree  well 
with  results  in  references  15  and  10  using  the  modified  mapping  collocation 
method.  Tracy,  in  a  private  communication,  pointed  out  a  discrepancy  of  five 
percent  for  the  case  of  four  cracks  with  crack  depth  c/t  =  0.5  in  a  tube 
having  outer  diameter  twice  that  of  inner  diameter. 

In  the  previous  work,l  the  finite  element  meshes  were  automatically 
generated  by  a  computer  program  for  various  values  of  N  and  c/t.  There  is  a 
possibility  of  excessive  distortion  of  elements  for  some  values  of  N  and  c/t. 
According  to  Sickles  and  Gif ford, 16  isoparametric  elements  suffer  a  loss  of 
accuracy  when  distorted  from  a  rectangular  shape.  They  recommended  a  '45 
degree  rule'  as  a  guide  to  construct  the  element  mesh.  Because  of  these 
reasons  and  our  interest  in  very  shallow  radial  cracks,  we  have  used  the 


Ipu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  For  a  Circular  Ring 
With  Uniform  Array  Radial  Cracks  Using  Isoparametric  Singular  Elements,” 
ASTM  STP-677 ,  1979,  pp.  685-699. 

l^Tracy,  P.  G.,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer 
and  Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol. 
11,  1979,  pp.  291-300. 

14Pu,  S.  L.  and  Hussain,  M.  A.,  "The  Collapsed  Cubic  Isoparametric  Element  as 
a  Singular  Element  for  Crack  Problems,"  International  Journal  of  Numerical 
Methods  in  Engineering,  Vol.  12,  1978,  pp.  1727-1742. 

^Bowie,  0.  L.  and  Freese,  C.  E.,  "Elastic  Analysis  For  a  Radial  Crack  in  a 
Circular  Ring,"  Journal  of  Engineering  Mechanics,  Vol.  4,  1972,  pp.  315-321. 

l^sickles,  J.  B.  and  Gifford,  L.  N.,  "A  Further  Study  of  Accuracy  Loss  in 
Distorted  Isoparametric  Finite  Elements,"  DTNSRDC  Report  M-50,  1979. 
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enriched  quadrilateral  elements^  as  crack  tip  elements  in  this  study,  A  typ¬ 
ical  finite  element  idealization  is  shown  in  Figure  1(a).  For  shallow  cracks, 
the  section  containing  the  crack  tip  changes  slightly  as  shown  in  Figure  1(b) . 
When  thermal  simulation  is  used  with  the  finite  element  method,  the  circle  r  = 
p  must  be  a  side  of  quadrilateral  elements  in  the  finite  element  idealization 
since  thermal  loads  are  different  in  the  two  regions  of  r  <  p  and  r  >  p.  The 
actual  finite  element  idealization  for  a  particular  geometry  and  thermal  loads 
may  be  mbdified  slightly  from  those  shown  in  Figure  1 . 

The  finite  element  computer  program  APES,  an  acronym  for  Axisymmetric/ 
Planar  Elastic  Structures,  is  used  for  all  finite  element  computations  in  this 
report.  This  powerful  computer  program  has  been  continuously  improved  with 
new  features  including  the  addition  of  thermal  loading  for  fracture 
analysis, The  APES  results  for  stress  intensity  factors  for  ID  radial 
cracks  for  a  cylinder  of  b  =  2  are  given  in  Table  I  for  three  types  of 
loading,  namely  (a)  uniform  tension  pQ  on  OD,  (b)  uniform  pressure  p^  on  ID 
(no  crack  face  pressure)  and  (c)  thermal  loading  equivalent  to  a  100  percent 

overstrain  residual  stress. 2  The  new  results  for  uniform  tension  on  OD  serve 

j 

as  a  check  to  the  previously  reported  results  using  collapsed  singular  crack 


^Hussain,  M.  A.,  Pu.  S.  L.,  Vasilakis,  J.  D.,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 

^Gifford,  L.  N.,  Jr.,  "APES  -  Second  Generation  Two-Dimensional  Fracture 
Mechanics  and  Stress  Analysis  by  Finite  Elements,"  DRNSRDC  Report  4799. 
1975. 

^Gifford,  L.  N. ,  Jr.,  "APES  -  Finite  Element  Fracture  Mechanics  Analysis: 
Revised  Documentation,"  DTNSRDC  Report  79/023,  1979. 
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tip  elements.1  The  new  results  agree  within  three  percent  with  those  reported 
in  references  15  and  10.  For  instance,  the  result  for  b=2,N=4,  c/t= 

0.5,  which  was  K(p0)/p0/iTc  =  2.990  in  reference  1,  is  now  K(p0)/p0/ttc  = 

3.149  which  agrees  with  Tracy's  Hj  =  1.181®  within  0.1  percent. 

An  explanation  is  in  order  regarding  negative  stress  intensity  factors 
shown  in  Table  I.  A  crack  remains  closed  in  a  compressive  residual  stress 
region.  The  stress  intensity  factor  is  zero.  The  crack  will  open  when  a  suf¬ 
ficiently  large  internal  pressure  is  applied.  The  negative  value  of  SIF  is 
convenient  in  measuring  the  crack  resistance  against  the  opening  by  internal 
pressure  and  it  should  be  understood  as  such. 

The  combination  of  the  finite  element  method  and  thermal  simulation 
method  can  be  used  to  compute  SIF  for  any  degree  of  partial  autofrettage . 
However,  it  is  expensive  and  time  consuming  to  use  finite  element  for 
parametric  studies. 

Therefore,  we  seek  alternative  methods  for  the  computation  of  SIF  for 
radial  cracks  in  a  partially  autofrettaged  tube  in  the  following  two  sections. 


iPu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  For  a  Circular  Ring 
With  Uniform  Array  Radial  Cracks  Using  Isoparametric  Singular  Elements," 
ASTM  STP-677 ,  1979,  pp.  685-699. 

^Tracy,  P.  G.,  "Elastic  Analysis  of  Radial  Cracks  Emanating  From  the  Outer 
and  Inner  Surfaces  of  a  Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol. 
11,  1979,  pp.  291-300. 

15Bowie,  0.  L.  and  Freese,  C.  E.,  "Elastic  Analysis  For  a  Radial  Crack  in  a 
Circular  Ring,"  Journal  of  Engineering  Mechanics,  Vol.  4,  1972,  pp.  315-321. 
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TABLE  I.  DIMENSIONLESS  SIF,  K(p0) /p0/irc,  K(  p^ /pi/rTc  and  Kc(  e=l)/o0/iTC 


OBTAINED  FROM 

APES  FOR 

A  CYLINDER  b  =  2 

FOR  VARIOUS 

N  AND  c/t 

N 

Loading 

c/ 1  =  0 . 

,05  c/t  =  0.1 

c/t  =  0.2 

c/t  =  0.3 

Po 

2.874 

2.825 

2.828 

2.890 

1 

Pi 

1.703 

1.711 

1.667 

1.646 

100%  OS 

-0.967 

-0.896 

-0.758 

-0.650 

Po 

2.891 

2.874 

3.014 

3.279 

2 

Pi 

1.762 

1.745 

1.766 

1.872 

100%  OS 

-0.980 

-0.919 

-0.813 

-0.745 

Po 

2.878 

2.860 

2.882 

3.034 

3 

Pi 

1.762 

1.731 

1.689 

1.728 

100%  OS 

-0.980 

-0.907 

-0.776 

-0.683 

Po 

2.866 

2.826 

2.782 

2.833 

4 

Pi 

1.764 

1.710 

1.629 

1.611 

100%  OS 

-0.981 

-0.895 

-0.745 

-0.628 

Po 

2.843 

2.753 

2.578 

2.504 

6 

Pi 

1.753 

1.665 

1.507 

1.418 

100%  OS 

-0.975 

-0.871 

-0.684 

-0.539 

Po 

2.792 

2.590 

2.217 

2.056 

10 

Pi 

1.721 

1.566 

1.291 

1.155 

100%  OS 

-0.957 

-0.817 

-0.575 

-0.417 

Po 

2.548 

2.062 

1.635 

1.538 

20 

Pi 

1.570 

1.243 

0.944 

0.854 

100%  OS 

-0.871 

-0.641 

-0.401 

-0.284 

Po 

2.224 

1.687 

1.363 

1.282 

30 

Pi 

1.367 

1.014 

0.782 

0.709 

100%  OS 

-0.756 

-0.515 

-0.322 

-0.228 

Po 

1.959 

1.465 

1.169 

1.124 

40 

Pi 

1.195 

0.879 

0.670 

0.621 

100%  OS 

-0.659 

-0.442 

-0.277 

-0.197 
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METHOD  OF  LOAD  RELIEF  FACTOR 


Baratta^  extended  the  coefficient  of  load  relief  of  Neuber^  to  estimate 
SIF  arising  from  multiple  cracking  in  a  thick-wall  cylinder.  He  defined  the 
ratio 

Km 

R  =  —  (12) 

Kn 

as  the  load  relief  factors,  with  Km  and  Kn  based  on  a  solution  due  to  Tweed 
and  Rooke.^O  Assuming  that  this  R  value  remains  nearly  constant  from  an 
infinitely  thick  cylinder  to  a  cylinder  with  finite  thickness,  he  estimated 
Km,f,  since  Kn>f  was  known  for  n  =  1  and  n  =  2  from  Bowie  and  Freese.^  His 
estimates,  based  on  the  crude  assumption,  have  a  discrepancy  as  high  as  20 
percent  with  our  finite  element  results  for  the  same  geometry  and  loading. 
Parker  and  Farrow^l  obtained  R  values  for  various  numbers  of  cracks  using  our 
results  in  reference  1.  They  assumed  that  R  values  vary  with  geometrical 
configurations,  but  are  independent  of  loading.  Therefore,  they  obtained  SIF 
for  N  ID  radial  cracks  for  a  100  percent  overstrained  tube  by  applying  the 


1Pu,  S.  L.  and  Hussain,  M.  A.,  "Stress  Intensity  Factors  For  a  Circular  Ring 
With  Uniform  Array  Radial  Cracks  Using  Isoparametric  Singular  Elements," 
ASTM  STP-677 ,  1979,  pp.  685-699. 

^Baratta,  F.  I.,  "Stress  Intensity  Factors  For  Internal  Multiple  Cracks  in 
Thick-Walled  Cylinders  Stressed  by  Internal  Pressure  Using  Load  Relief 
Factors,"  Engineering  Fracture  Mechanics,  Vol.  10,  1978,  pp.  691-697. 
l^Bowie,  0.  L.  and  Freese,  C.  E.,  "Elastic  Analysis  For  a  Radial  Crack  in  a 
Circular  Ring,"  Journal  of  Engineering  Mechanics,  Vol.  4,  1972,  pp.  315-321. 
l^Neuber,  h.,  "Theory  of  Notch  Stresses,"  AEC  TR  4547,  1958. 

^Tweed,  J.  and  Rooke,  D.  P.,  "The  Stress  Intensity  Factor  For  a  Crack  in 
Symmetric  Array  Originating  From  a  Circular  Hole  in  an  Infinite  Solid," 
Journal  of  Engineering  Science,  Vol.  13,  1975,  pp.  653-662. 

-^Parker,  A.  P.  and  Farrow,  J.  R. ,  "Stress  Intensity  Factors  For  Multiple 
Radial  Cracks  Emanating  From  the  Bore  of  an  Autofrettaged  or  Thermally 
Stressed  Thick  Cylinder,"  Materials  Branch  Technical  Note  MAT/ 20,  Royal 
Military  College  of  Science,  England,  1979. 
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load  relief  factors  to  the  solution  of  SIF  for  a  100  percent  autof rettaged 
tube  with  two  ID  radial  cracks  due  to  Grandt.^2 

If  the  load  relief  factor  works  for  a  fully  autofrettaged  tube,  then  it 
should  work  for  a  partially  overstrained  tube  also.  Our  study  reveals  that  R 
values  do  vary  with  the  nature  of  load.  Strictly  speaking,  the  concept  of 
load  relief  factor  does  not  work.  However,  for  a  thick  wall  cylinder  subjec¬ 
ted  to  the  following  three  types  of  loading:  (a)  uniform  tension  p0  on  OD, 

(b)  uniform  pressure  p^  on  ID,  and  (c)  residual  stress  due  to  100  percent  over¬ 
strain,  R  values  vary  within  ±  five  percent  except  for  large  values  of  N  and 
c/t.  Table  II  gives  R  values  for  various  values  of  N  and  c/t  under  these 
three  types  of  loading.  From  Table  II  it  can  be  seen  that  R  values  remain 
nearly  constant  for  any  given  N  and  c/ t  no  matter  if  the  loading  is  pQ  on  OD  or 
p^  on  ID.  This  property  enables  us  to  estimate  fairly  accurately  the  SIF  for 
Pi  from  values  for  pQ  and  vice  versa.  The  R  values  corresponding  to  the  100 
percent  overstrain  residual  stress  agree  within  two  percent  with  R  values  for 
p0  if  N  <  4  and  c/ t  <  0.3.  Therefore,  SIF  may  be  estimated  for  N  <  4,  c/t  < 

0.3  for  partially  autofrettaged  cases  from  R  values  for  pQ  and  from  known 
results  of  SIF  for  N  =  2  for  the  same  crack  depth  c/t  and  loading.  The  error 
may  exceed  ten  percent  if  the  load  relief  factor  method  is  used  to  estimate 
SIF  for  cracked  tubes  with  autofrettage  residual  stress  from  R  values  obtained 
for  pQ  on  OD  for  N  >  10  and  c/t  >  0.2.  For  loading  conditions  other  than  the 
three  types  mentioned,  the  method  of  load  relief  factor  should  be  applied  with 
care. 

^^Grandt,  A.  F.,  "Two  Dimensional  Stress  Intensity  Factor  Solutions  For 

Radially  Cracked  Rings,"  Technical  Report  AFML-TR-7 5-121 ,  Air  Force  Materials 
Laboratory,  Wright-Patterson  Air  Force  Base,  1975. 
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TABLE  II.  LOAD  RELIEF  FACTOR  R  =  Kn/Kn=2  FOR  A  CYLINDER  OF  b  =  2 
SUBJECTED  TO  THREE  TYPES  OF  LOADING 


N 

Loading 

c/t  =  0 .05 

c/t  =  0.1 

c/t  =  0.2 

c/t  =  0.3 

Po 

0.994 

0.983 

0.938 

0.881 

1 

Pi 

0.967 

0.980 

0.943 

0.897 

100%  OS 

0.987 

0.976 

0.932 

0.873 

Po 

0.996 

0.995 

0.956 

0.925 

3 

Pi 

1.000 

0.992 

0.956 

0.923 

100%  OS 

1.000 

0.987 

0.954 

0.918 

Po 

0.992 

0.983 

0.923 

0.864 

4 

Pi 

1.002 

0.980 

0.922 

0.860 

100%  OS 

1.002 

0.974 

0.917 

0.844 

Po 

0.984 

0.958 

0.855 

0.764 

6 

Pi 

0.995 

0.954 

0.853 

0.757 

100%  OS 

0.995 

0.948 

0.841 

0.723 

Po 

0.966 

0.901 

0.735 

0.627 

10 

Pi 

0.977 

0.897 

0.731 

0.617 

100%  OS 

0.977 

0.889 

0.707 

0.560 

Po 

0.882 

0.718 

0.543 

0.469 

20 

Pi 

0.891 

0.713 

0.534 

0.456 

100%  OS 

0.889 

0.697 

0.493 

0.381 

Po 

0.769 

0.587 

0.452 

0.391 

30 

Pi 

0.776 

0.581 

0.443 

0.379 

100%  OS 

0.772 

0.561 

0.396 

0.307 

Po 

0.678 

0.510 

0.388 

0.343 

40 

Pi 

0.678 

0.504 

0.379 

0.331 

100%  OS 

0.672 

0.481 

0.340 

0.264 
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WEIGHT  FUNCTION  METHOD 


Bueckner^  and  Rice^  have  shown  that  knowledge  of  the  SIF  and  dis¬ 
placement  field  for  a  flaw  geometry  enables  construction  of  a  weight  function 
which  depends  only  on  geometry.  With  the  weight  function  method,  one  may 
obtain  SIF  for  any  other  symmetric  loading  applied  to  the  same  geometry. 

Grandt  has  applied  this  technique  to  obtain  SIF  for  a  large  plate  containing 
radial  hole  cracks?  and  for  radially  cracked  rings22  loaded  with  arbitrary 
symmetric  crack  pressure.  The  SIF  K  for  a  specified  crack  face  loading  pc(x), 
based  on  the  weight  function  approach,  is  given  by 

H  .c  3v 

K  »  /  pc(x)  --  dx  (13) 

K*  o  3c 

where  H  is  a  constant,  H  =  E  for  plane  stress  and  H  =  E/(l-v2)  for  plane 
strain,  K*  is  the  known  SIF  for  a  given  loading  applied  to  the  flaw  geometry 
of  interest,  x  is  the  distance  from  the  edge  of  the  hole,  v  is  the  crack 
opening  profile  corresponding  to  the  known  SIF  K* .  In  Eq.  (13)  the  only 
undefined  term  is  the  partial  derivative  3v/3c.  The  finite  element  results  of 
the  y-component  of  displacement  at  nodal  points  along  the  crack  face  are  not 
enough  for  the  determination  of  the  crack  profile.  Grandt  used  the  assumption 


^Bueckner,  H.  F.,  "A  Novel  Principle  For  the  Computation  of  Stress  Intensity 
Factors,"  Z.  Agnew.  Math.  Mech.,  Vol.  50,  1970,  pp.  529-546. 

^Rice,  J.  R. ,  "Some  Remarks  on  Elastic  Crack-Tip  Stress  Fields,"  Int.  Journal 
of  Solids  and  Structures,  Vol.  8,  1972,  pp.  751-758. 

? Grandt,  A.  F.,  "Stress  Intensity  Factors  For  Some  Through-Cracked  Fastener 
Holes,"  Int.  Journal  of  Fracture,  Vol.  11,  1975,  pp.  283-294. 

22Grandt,  A.  F.,  "Two  Dimensional  Stress  Intensity  Factor  Solutions  For 
Radially  Cracked  Rings,"  Technical  Report  AFML-TR-7 5-121 ,  Air  Force 
Materials  Laboratory,  Wright-Patterson  Air  Force  Base,  1975. 
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of  conic  sections  due  to  Orange. 6  Andrasic  and  Parker^  used  the  method  of 
virtual  crack  extension* 2^  and  B-spline  curve  fitting.  In  this  report  the 
loading  of  the  cracked  tube  is  limited  to  a  combination  of  internal  pressure 
and  the  autofrettage  residual  stresses.  A  set  of  algebraic  equations  is  used 
in  lieu  of  the  determination  of  9v/9c. 

The  crack  face  loading  pc(x)  in  Eq.  (13)  for  a  cracked  tube  subjected  to 
the  autofrettage  residual  stress  is  given  by  the  hoop  stress,  Eqs.  (3)  and 
(4),  for  an  uncracked,  overstrained  tube 


1 


pc(x)  oe(x) 


-< 


--  [(2-Pi)  -  PxCl+x)-2  +  2  log  ( 1+x) ] 
/3 


—  [(p2-Pl)b  2  +  (p2-Pi)(l+x)  2] 
/3 


0  <  x  <  et  (14) 


et  <  x  <  t  (15) 


For  e  =  1,  substituting  from  Eq.  (14)  into  Eq.  (13),  the  following  is 


obtained : 


Kc(e=l)  ! 

-  =  -  [{2  -  PL(b)}  Kc(l)  -  Px(b)  Kc( r-2)  +  2KC( log  r) ]  (16) 

oo/ttc  /3ttc 


where 

H  c  9v 

Kc(l)  *  —  /  —  dx  (17) 

K*  o  9c 


^Orange,  T.  W.,  "Crack  Shapes  and  Stress  Intensity  Factors  for  Edge-Cracked 
Specimens,"  ASTM  STP-513,  1972,  pp.  71-78. 

^^Andrasic,  C.  P.  and  Parker,  A.  P.,  "Weight  Functions  For  Cracked  Curved 
Beams,"  Second  International  Conference  on  Numerical  Methods  in  Fracture 
Mechanics,  Swansea,  U.K.,  1980. 

2^Parks,  D.  M.  and  Kamenetzky,  E.  M.,  "Weight  Functions  From  Virtual  Crack 
Extensions,"  International  Journal  For  Numerical  Methods  in  Engineering, 
Vol .  14,  1979,  pp.  1693-1706. 
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(18) 


9  H  .c  9  3v 

Kc(r“2)  =  —  /  (1+x)"2  —  dx 
K*  o  3c 


H  ,c  3v 

Kc(log  r)  =  —  /  [ log( 1+x) ]  --  dx 

K*  o  3  c 


(19) 


These  may  be  termed  as  functional  intensity  factors. 

From  the  Lame^  solution,  the  hoop  stress  in  an  uncracked  cylinder  subjected  to 
uniform  tension  p0  on  OD  is 

°e  b2  l 

(20) 


°e  b2  l 

—  =  (l  +  -^) 


Po  b2-l  r“ 

The  same  stress  under  internal  pressure  p-£  is 

°0  1  b2 

—  =  (1  +  -,)  (21) 

Pi  b2-l  r2 

Substituting  o@  for  pc  in  Eq.  (13),  the  SIF  for  a  radially  cracked  cylinder 
subjected  to  pQ  or  p^  is  given  by  one  of  the  following 


K(p0)  b2  KC(D  b2  Mr-2) 
Po/^  "  b2-l  /7Z  b2-l  /iH 


(22) 


K(pi)  1  Kc(l)  b2  Kc(r-2) 

Pi/7H  ~  b2-l  b2-l 


(23) 


If  the  left  hand  sides  of  Eqs.  (22)  and  (23)  are  numerically  determined  from 
the  finite  element  method,  then  Kc(l)  and  Kc(r-2)  can  be  solved.  Inserting 
these  values  into  Eq.  (16),  Kc(log  r)  can  be  solved  if  Kc(e=l)  is  known  from 
the  finite  element  computation. 

Once  Kc(l),  Kc(r~2) ,  and  Kc(log  r)  are  known  for  a  given  geometry  (b,  N, 
c/ t  are  fixed),  the  SIF  for  the  flawed  cylinder  subjected  to  a  residual  stress 
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corresponding  to  a  given  e  (*1)  can  be  computed  from 
Kc(e)  i 

"  — 1  [  {2-P1(p)}Kc(l)  -  P1(p)Kc(r"2)  +  2KC( log  r)  ]  (24) 

a0/7Tc  /3ttc 

provided  c  <  et. 


In  case  c  >  et,  the  crack  face  pressure  is  a  combination  of  Eqs.  (14)  and 
(15).,  therefore,  Eq.  (24)  is  not  valid.  If  c/ 1  =  e+5  and  0(6  ((  1,  we  may 
use  Eq.  (24)  to  compute  an  approximate  SIF  and  then  use  the  following  equation 
to  compute  a  corrective  SIF,  Kg 


where 


1  H 
/3tFc  K* 


/ 


( e+6 ) t 


et 


Pc(x) 


9v 

—  dx 
9c 


(25) 


Pc(x)  “  (“1+2  log  p)  +  p2(l+x)  2  -  2  log(l+x)  (26) 
The  final  result  of  SIF  for  a  small  6  >  0  is  the  sum  of  Kc(e),  Eq.  (24),  and 
Kfi,  Eq.  (25). 


An  approximate  method  for  Kg  is  based  on  the  Westergaard  near  field 
solution.25  The  crack  opening  displacement  v( 5)  near  a  crack  tip  due  to  an 
arbitrary  load  is  given  in  terms  of  SIF  K*  at  the  crack  tip  by 

1/2 


2K*  25 

v(5)  = - (“) 

H  tr 


The  derivative  with  respect  to  the  crack  length  is 

9  v  K*  2  1//2  ,  . 

=  —  (-)  (r1/2  +  51/2/c) 

9c  H  t: 


(27) 


(28) 


)  c 

Westergaard,  H.  M.,  Bearing  Pressures  and  Cracks,"  Transactions  of  the 
ASME,  Journal  of  Applied  Mechanics,  1939. 
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The  length  variable  £  is  defined  by 

5  =  -(x-c)  (29) 

Substituting  from  Eqs.  (26)  and  (28)  into  (25)  and  using  Eq.  (29),  we  have 
KS  l 

-  =  -  /2/tT  {(_!  +2  iog  p)(ll  +  Xl»)  +  p2(x2  +  l2«)  _  2(I3  +  I3')}  (30) 

o0/ttc  /3ttc 


where 


IX  =  2 /fit  ,  Ii'  =  —  ( 6t)3/2 

3c 

1  (i+c)_1/2 

I2  =  7“  t - log  D(  p)  ] 

1+c  p  2 

1  St*  (l+c)_1/2 

i2l  =  -  [ - + - log  D(  p) ] 

c  p  2 


(31) 


(32) 


/  1+c  -  /fit 

1 3  =  -2  [  ( 2-log  p)/fit  +  (l+c)l/2  log - ] 

/  1+c  +  /fic 


I3'  =  —  [ /fit( 6 1{ log  p  -  2( 1+c)  -  2fit/3) }  -  (1+c)3/2  log  D( p) ] 
3c 


(33) 


with 


D(p)  =  [2(l+c)  -  p  -  2/ fit( 1+c) ] / p 


(34) 


It  should  be  noted  that  Eq.  (30)  is  independent  of  N.  It  works  for  small  N 
and  5.  If  N  is  small,  it  may  work  for  a  relatively  larger  6.  But  when  N  is 
large,  the  crack  interaction  is  strong,  6  must  be  small. 
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NUMERICAL  RESULTS 


Using  Eqs.  (16),  (22),  and  (23)  and  finite  element  results  in  Table  I,  we 
obtain  values  of  Kc(p)/p/irc,  Kc(pr“2)/p/Trc,  and  Kc(p  log  r)/p/lTc.  Figures  2 
through  4  are  plots  of  these  values  as  a  function  of  c/ t  for  various  values  of 

N.  The  use  of  these  graphs  and  Eq.  (24)  gives  Kc(e)/a0/Trc  for  any  value  of  e, 
et  >c,  for  a  given  geometry.  A  graph  of  Kc(  e)/a0/rirc  versus  e  is  shown  in 
Figure  5  for  various  c/t  and  for  N  =  2  and  N  =  40.  Figure  6  is  another  way  of 
presenting  Kc(e)/a0/7rc  in  which  e  is  fixed  but  N  varies. 

The  numerical  results  given  previously  are  enough  for  an  estimate  of  SIF 
for  any  assigned  values  of  N,  c/t,  and  e.  For  example,  if  the  SIF  is  desired 
for  N  =  8,  c/t  -  0.15,  and  e  =  0.75,  readings  are  taken  from  Figures  2  through 
4  as  follows:  Kc(p)/p/ttc  =  1.03,  Kc(pr-2)/p/iTC  =  0.86,  and  Kc(p  log  r)/p/7rc  = 

O. 091.  For  e  =  0.75  we  have  from  Eqs.  (5)  and  (24)  Kc(  e  =  0.75)/oo/ttc  = 

-0.67.  In  another  example,  if  the  SIF  for  N  =  2 ,  c/ 1  =  0.30  and  e  =  0.25  is 
desired,  we  first  compute  Kc(e  =  0.25)/ao/TTC  from  Eq.  (24)  with  Kc(1)//ttc  - 
1.41  from  Figure  2,  Kc(r“2)/7rc  =  1.05  from  Figure  3,  and  Kc(log  r)//iTc  =  0.22 
from  Figure  4.  The  result  from  Eq.  (24)  is  Kc(  e  =  0.25)/ao/7TC  «  -0.12.  In 
this  case  since  c/ 1  >  e ,  we  have  to  compute  K^/Oq/ttc  from  Eqs.  (30)  through 
(34)  with  6t  =  0.05.  The  corrective  SIF  is  K5/a0/7rc  -  -0.023  and  the  desired 
SIF  is  Kc(e  =  0.25)/oo/7tc  =  -0.143.  A  finite  element  computation  is  performed 
for  the  case;  the  result  is  also  -0.143.  If  in  the  previous  example  e  =  0.20 
is  desired,  i.e.,  <5t  =  0.1,  Eq.  (24)  gives  Kc/ o0/ttc  =  -0.0209  and  Eq.  (30) 
gives  K5/o0/ttc  =  -0.0668.  The  final  result  is  Kc(  e  =  0.2)/oo/7tc  =  -0.088 
which  is  5.6  percent  less  than  the  finite  element  result  Kc(e  =  0.2)/oo/ttc  = 

-0.093.  This  indicates  that  6  must  be  fairly  small  for  Eq.  (30)  to  be  valid. 
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For  the  geometry  N  =  10,  c/t  =  0.3,  the  finite  element  computation  gives  Kc(e 
=  0.3)/ao/iTC  =  -0.079  and  Kc(e  =  0.25)/oo/ttc  =  -0.044.  The  corresponding 
values  computed  from  Eq.  (24)  and  from  Eqs.  (24)  and  (30)  are  -0.080  and 
-0.0477. 

The  stress  intensity  factors  due  to  a  combination  of  residual  stress  and 
internal  pressure  on  ID  and  on  crack  faces  can  be  readily  computed  by  linear 
superposition.  If  the  applied  pressure  is  p^  =  aQ/ f ,  where  f  is  a  constant, 
the  resultant  SIF’s  for  f  =  1.5  and  3,  for  N  =  1  and  40,  and  for  e  =  1.0  are 
shown  in  Figure  7.  We  purposely  keep  the  negative  SIF  for  shallow  cracks  In 
the  case  of  f  =  3  and  e  =  1.0.  The  correct  SIF  should  be  zero,  which  means 
that  a  single  crack  or  a  set  of  multiple  cracks  remains  closed  due  to  the  high 
compressive  residual  stress  near  the  bore.  The  negative  SIF  gives  a  little 
more  information  than  SIF  =  0.  The  graph  also  shows  a  higher  SIF  for  N  =  1 
than  that  for  N  =  40  for  a  given  e  and  f.  One  may  conclude  that  it  is  worse 
to  have  a  single  radial  crack  than  to  have  a  large  number  of  radial  cracks  in 
an  autof rettaged  tube.  This  conclusion  is  similar  to  that  found  in  the  study 
of  non-autof rettaged  tubes.  With  some  slight  modification,  the  method  used 
here  can  be  applied  to  multiple  OD  cracks.  The  modifications  and  numerical 
results  for  OD  cracks  are  to  be  reported  elsewhere. 
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CONCLUSIONS 


The  finite  element  method  together  with  the  thermal  simulation  method  can 
be  used  to  compute  the  stress  intensity  factor  for  multiple  radial  cracks 
emanating  from  the  bore  of  a  partially  autofrettaged  tube.  The  finite  element 
results  of  K(p0)  ,  K(p^),  and  Kc(e  =  1.0)  can  lead  to  a  system  of  algebraic 
equations  for  solving  Kc(l),  Kc(r“2),  and  Kc(log  r) .  Using  these  results  the 
weight  function  concept  gives  an  alternative  method  for  the  determination  of 
,SIF  for  any  degree  of  partial  autofrettage  provided  that  the  crack  depth  c/t 
is  not  greater  than  the  percentage  of  autofrettage  e.  A  correction  formula  is 
supplied  for  the  SIF  when  c/t  is  slightly  greater  than  e  when  N  is  small. 

These  expressions  yield  quite  accurate  results  and  can  save  a  great  deal  of 
computing  time. 

The  useful  tube  life  is  prolonged  because  of  the  autofrettage  process. 

The  cylinder  with  two  diametrically  opposed  cracks  remains  in  general  the 
weakest  configuration.  For  more  than  two  cracks,  the  stress  intensity  factor 
decreases  as  the  number  of  cracks  increases.  The  fatigue  life  is  longer  for 
cylinders  with  more  radial  cracks. 
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Figure  1(a). 


(b) 


Figure  1(b). 


CRACK  TIP 


A  typical  finite  element  idealization. 


CRACK  TIP 


Idealization  for  very  shallow  cracks. 
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P\fnc 


Figure  2.  Stress  intensity  factors  as  a  function  of  c/t  for  N  radial 
cracks  subjected  to  the  crack  face  pressure  pc(x)  =  p. 
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Kc(pr'2) 


Figure  3.  Stress  intensity  factors  as  a  function  of  c/t  for  N  radial 

cracks  subjected  to  the  crack  face  pressure  pc(x)  =  p(l+x)~2. 
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Figure  4. 


Stress  intensity  factors  as  a  function  of  c/t  for  N  radial 
cracks  subjected  to  the  crack  face  pressure  pc(x)  =  p  log  (1+x) . 
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Figure  5.  Stress  intensity  factors  as  a  function  of  e  in 
an  autofrettaged  cylinder  of  b  =  2. 
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Figure  6.  Stress  intensity  factors  as  a  function  of  N  in 
a  fully  autofrettaged  cylinder  of  b  =  2. 


Figure  7.  Stress  intensity  factors  in  an  autofrettaged  cylinder  of  b  =  2 
subjected  to  internal  pressure  cr0/f  on  ID  and  on  crack  faces. 
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